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Abstract 

Valence u- and d-quarks distributions in proton are calculated in QCD at low and 
intermediate x, basing on the operator product expansion (OPE). The imaginary part of 
the virtual photon scattering amplitude on quark current with proton quantum numbers 
is considered. The initial and final virtualities and of the currents are assumed 
to be large, negative and different, 7^ pI- OPE in P2 to dimension 6 

operators was performed. Double dispersion representations inpf of the amplitudes 
in terms of physical states contributions are used. Putting them to be equal to those 
calculated in QCD, the sum rules for quark distributions are found. The double Borel 
transformations are applied to the sum rules. Leading order perturbative corrections 
are accounted. Valence quark distributions are found: u(x\ at 0.15 < x < 65, fi(x)„ 
at 0.25 < X < 0.55 with an accuracy ~ 30% in the middles and ~ 50% at the ends of 
these intervals. The quark distributions obtained are in agreement with those found 
from the analysis of hard processes data. 



FACS: ll.15.Tk; 11.55.Hx; 12.38.Lg. 
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1. Introduction 



Quark and gluon distributions in hadrons are not fully understood in QCD. QCD predicts 
the evolution of these distributions with in accord with the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) [l]-[3] equations, but not the initial values from which this evolution 
starts. The standard way of determination of quark and gluon distributions in nucleon is 
the following [4-8] (for the recent review see [9]). At some = Ql (usually, at low or 
intermediate ~ 2 — bGeV"^) the form of quark (valence and sea) and gluon distributions 
is assumed and characterized by the number of free parameters. Then, by using DGLAP 
equations, quark and gluon distributions are calculated at all and x and compared with 
the whole set of the data on deep inelastic lepton-nucleon scattering (sometimes also with 
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prompt photon production, jets at high p± etc). The best fit for the parameters is found 
and, therefore, quark and gluon distributions are determined at all Q^, including their initial 
values q{Ql,x), g{QQ,x). Evidently, such an approach is not completely satisfactory from 
theoretical point of view - it would be desirable to determine the initial distribution directly 
from QCD. Also, if the form of initial distribution is not divined properly, then the fitting 
procedure will give the relative minimum, but not the absolute ones and the results, especially 
at low Ql could be wrong. (This danger does not, probably, exists for nucleon, but such a 
situation happens, e.g., for transversally polarized p-meson, where unusual form of initial 
valence quark distribution was found [10]). Finally, extrapolation from high Q^, to which the 
main part of the data belongs, to low is not a quite stable procedure and may introduce 
some errors. For all these reasons it is desirable to find quark and gluon distribution in 
hadrons at low ~ 2 — 5GeV'^ basing directly on QCD. 

In this paper we calculate u and d valence quark distributions in proton. The idea of 
the method was suggested in [11] and developed in [12-14]. Recently, the method had been 
improved and valence quark distributions in pion [15] and transversally and longitudinally 
polarized p-meson [10] had been calculated, what was impossible in the initial version of the 
method. The idea of the approach (in the improved version) is to consider the imaginary 
part (in s-channel) of a four-point correlator n(pi,p2, Q, q') corresponding to the non-forward 
scattering of two quark currents, one of which has the quantum numbers of hadron of interest 
(in our case - of proton) and the other is electromagnetic (or weak). It is supposed that 
virtualities of the photon q'^,q'^ and hadron currents Pi,P2 are large and negative = 
^ IpII ^ ^c^j where Rc is the confinement radius. It was shown in [12] that in this 
case the imaginary part in s-channel [s — {pi + qY] of n(pi,p2; qi, q') is dominated by a small 
distance contribution at intermediate x. (The standard notation is used: x is the Bjorken 
scaling variable, x = — g^/2z/ , u = piq). The proof of this statement is given in ref. [12]. So, 
in the mentioned above domain of g^, g'^, Pi,P2 and intermediate x Imll{p,p2; q, q') can be 
calculated using the perturbation theory and the operator product expansion in both sets 
of variables q^ — g'^ and pi, pi- The approach is inapplicable at small x and x close to 1. 
This can be easily understood for physical reasons. In deep inelastic scattering at large \q^\ 
the main interaction region in space-time is the light-cone domain and longitudinal distances 
along the light-cone arc proportional to 1/x and become large at small x [16, 17]. For OPE 
validity it is necessary for these longitudinal distances along light-cone to be also small, that 
is not the case at small x. At 1 — x <C 1 another condition of applicability of the method is 
violated. The total energy square s = Q'^{l/x — 1)+Pi = —q^ is not large at 1 — a; -C 1. 
Numerically, the typical values to be used below are ~ SGeV^^, p\ ~ —IGeV^. Then, 
even at a; ^ 0.7, s ~ IGeV^, i.e., at such x we are in the resonance, but not in the scaling 
region. So, one may expect beforehand, that our method could work only up to a; ~ 0.7. The 
inapplicability of the method at small and large x manifests itself in the blow-up of higher 
order terms of OPE. More precise limits on the applicability domain in x will be found from 
the magnitude of these terms. 

The further procedure is common for QCD sum rules. On one hand the four-point cor- 
relator Il{pi,p2;q,q') is calculated by perturbation theory and OPE. On the other hand, the 
double dispersion representation in pl,pl in terms of physical states contributions is written 
for the same correlator and the contribution of the lowest state is extracted using the Borel 
transformaion. By equalling these two expression the desired quark distribution is found. 

In Sec. 2 the general outline of the method is presented. It is explained here, why it is 
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necessary to consider at the beginning the nonequal pi and pi and only at the end of the 
calculation to go to the forward scattering amphtude. Sec. 3 presents the results of the bare 
loop calculation (d — term in OPE) and the LO perturbative corrections to it. Sec.4 
is devoted to the calculation of power corrections to the QCD side of the sum rules: the 
contribution of gluonic condensate (d=4), as{il)ip)^ term (d=6). The valence u and ri-quark 
distributions are numerically calculated in Sec. 5, the errors are estimated and the results are 
compared with those obtained from the deep inelastic scattering data using the evolution 
equations. Sec. 6 contains our conslusion. 



2. The outline of the method 

Consider the 4-current correlator which corresponds to the virtual photon scattering on 
the quark current with quantum number of proton: 

T'"'{pi,P2,q,q') = -iJ d^xd^yd^z ■ e*(pi^+w-P2.) 

•(0|T{?7(x), r/{ylf/{Q\ (1) 

where r^(a;) is the three-quark current (loffe current [18]). Choose the currents in the form 
— ^liJ.'^-i — d'y^d, i.e. as an electromagnetic current which interacts only with u{d) 
quark (with unit charges). Such a choice allows us to get sum rules separately for distribuion 
functions of u and d quarks. The general method of calculation of distribution functions 
from consideration of the 4-point correlator in the QCD sum rules for light quarks was used 
in [12]. Unfortunately, this approach being in the form used in [12] did not give completely 
satisfactory sum rules for quark distribution functions (especially for d-quark in nucleon) and 
did not allow one to calculate quark distributions in mesons. Recently, in ref. [10,15] was 
suggested a generlized method of calculation of quark distribution functions, which made it 
possible to get much more reliable sum rules, what, in particular, made it possible to find 
valence quark distributions in tt and p- mesons. The main difference of this method is that the 
hadronic current momenta are put to be unequal, pf ^ pi and then independent borelization 
over p1 and p^ is performed and only at the very end the Borel the parameters Mf and M| 
are put to be equal. (Unlike the approach of ref. [12] where pi was put to be equal to p2 and 
single borelization is performed from the very beginning). The described procedure allows 
one to kill nondiagonal transitions of the type 

(o|/|/^*)(/^*|j?(y)j:'(o)|/i)(/i|/|o) (2) 

and thus makes it possible to separate the diagonal transition of interest 

(0|j''|/^)(/^|^(y)j:'(0)|/^)(/.|j'lO), (3) 

without using additional tricks like differentiation over Borel parameter M^, what strongly 
worsens the accuracy of the sum rules. 

Let us briefly remind the main points of derivation of the sum rules in the method under 
investigaion (for details see [10, 15]. The most general form of the double dispersion relation 
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(in pf and P2) for imaginary part of the correlator (1) (we omit for some time all indeces 
assuming the necessary invariant amplitude to be chosen) has the form 



ImT{pl,pl q\ s) = a{q\ s) + [ "^^"^^'f du + [ 

J u-pi J u-pi 



00 00 



+ / *M«,-P5)(».-Pi) 

(without loss of generality one may put = q'"^ , t = {pi — P2f' — The double Borel 
transformation in p\ and p| eliminates three first terms and we have 



00 00 
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BM^BMiImT{pl,pl,q^,s) = J dui J du2p{q'^,s,ui,U2)exp 



r2 ^-^A 



(5) 



where and are the squared Borel mass. The integration region with respect to ui, U2 
may be divided into four areas 

I. ui < So, U2 < So; 

II. ui < So, U2 > So; 

III. Ui > So, U2 < So; 

IV. Ml, U2 > Sq. 

Here sq is the continuum threshold in the standard QCD sum rule model of the hadronic 
spectrum with one lowest resonance plus continuum. Area I obviously corresponds to the 
resonance contribution and spectral density in this area can be written as 

p{ui,U2,x,Q'^) = gl ■ 27rF2{x,Q'^)6{ui - ml)6{u2 - ml), (6) 

where Qh are coupling constants of the corresponding hadronic current. As for the areas II- 
IV, they are exponentially suppressed, and, using the standard hypothesis of quark-hadron 
duality, we may estimate them as a bare loop contribution in the same integration region. 
Thus, the general form of the sum rule is 

ImT^CD + Power correction = 27rF2(x, Q^)gle " ^ ^ 

ImT^CD = J J p\u^,U2,x)du^du2e ^ (7) 



where one should go to the limit = M|. 

We choose — M| = 2M^ following [19] where it was shown that the value of the Borel 
mass square in the double sum rules is approximately twice as large than in single ones. 

Choose now invariant amplitude for the case of the proton current of interest in the 
correlator (1) (cf.[12]). Its imaginary (in s) part can be written as 
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Pi-m 

X Im{-i I d^xe-^'^y ■ {pi,r\T{j,iy),jM}\p2,r')} ■ v^' (p,)^^ (8) 

where v^{p) is proton spinor with polarization r and momentum p, X^is the couphng constant 
of proton with the current {0\ri\p,r) = Xnv'^{p). 

In order to choose the most suitable invariant amplitude, rewrite eq.(8) in the form 

.('r(p0/mf^j>t,''fc)")r(p2) (9) 

where Im Tj^ is the amplitude before averaging in proton spin, m is the proton mass. 
The general form of this amplitude is 

Im flP){p,,p2) = CiP^P. + C2(P^7. + ^/.7.) + C^PP^P, + ... 

+ (terms with r) (10) 

where P — 2i±£2 . f — — 

Let us now take into account that we are interested in such a combination of invariant 
amplitudes which in the spin-averaged matrix element 

$:r(pi)7mTj^>^(p2) (11) 

r 

appears at the kinematical structure P>^P'' (since this structure in the limit Pi — > P2 = P 
transforms into p^p'^, the coefficient at which is just ^2(0;)). 
Using the equation of motion 

Pi,2'"''{Pi,2) = \iPi,2 + m)v''{pi^2) (12) 

and 

VaiPl,2)vliPl,2) = (Pl,2 + m)a(3 (13) 

r 

it can be seen that the combination of invariant amplitudes in eq.(lO), which appears in eq.(9) 
at the kinematical structure PP'^P'^ coincides (up to numerical factor) with the combination 
of invariant functions at the structure P^Pu in (11). 

Thus, we come to a conclusion that the sum rules should be written for invariant amplitude 
at the kinematical structure pp/^p^ (in what follows we will denote it ImT/ v) . 

So, the sum rules for nucleon have the form 



27r A- 



2 



^ xg"''^(x)e-'"'/^' = ImT^^ + Power corrections (14) 



4M4 327r4 

where = 327r^A^; q^''^{x) are distribution functions of u{d) quark in nucleon, 7mT° is 
perturbative contribution, i.e. of a bare loop with perturbative corrections. (The continuum 
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contribution, i.e., of II, III, IV regions should be subtracted from IniT. Note, that really, 
the contribution of regions II and III to the bare loop is zero, since p° ~ 6{ui — U2)). 



3. Bare loop contribution and leading order perturbative corrections. 



Bare loop contribution to the sum rules is represented in Fig.l. In the calculation we use 
the technique described in ref.15. Appendix. The following formulae are exploited: 

du- 



{pi — ky{p2 — ky 8ra J {Pi — '^J') (pI ~ 

(The terms, vanishing at the double Borcl transformation are omitted). The results after the 
double Borel transformation arc the same as in the case for equal pi = p2 [12]: 

Im^ia) - E2{so/M') (16) 

where 

cp^{x) ^x{l-xy{l + 8x), ip^{x)^x{l-xf{l-2x), (17) 
So is the continuum threshold 

E^iz) ^l-{l + z + zy2)e-' (18) 
The substitution of eq.(16) into the sum rules (14) results in 

xq{xW = {x) ■ Mjp) (19) 

In the bare loop approximation the moments of the quark structure function are equal to 



1 

j q^{x)d 



'x = ^ Eo 







/Me ' 
q^{x)dx = 2 ^2 ^2 (20) 







Making use of relation \\e 1^ — M^E2 which follows from the sum rule for the nucleon 
mass (see [18]) in the same approximation, we get 
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1 




(21) 





In the bare loop approximation there also appears the sum rule for the second moment: 



Analogously to [12] one can show that relations (21), (22) hold also when taking into ac- 
count power corrections proportional to the quark condensate square in the sum rules for the 
4-point correlator (Fig.2) and in the sum rules for the nucleon mass. Relations (21) reflect 
the fact that proton has two w-quarks and one d-quark. Relation (22) expresses the momen- 
tum conservation law - in the bare loop approximation all momentum is carried by valence 
quarks. Therefore, the sum rules (21), (22) demonstrate that the zero order approximation is 
reasonable. In the real physical theory the regions x <^ 1 and 1 — x <^ 1 are off the frames 
of our consideration. However, in the noninteracting quark model which corresponds to the 
bare loop approximation, the whole region < x < 1 should be considered and relations 
(21), (22) should take place. 

Let us calculate the perturbative corrections to bare loop and restrict ourselves by the 
leading order (LO) corrections proportional to InQl/fj,^, where Ql is the point, where the 
quark distributions q{x, Qq) is calculated and //^ is the normalization point. In our case it is 
reasonable to choose /s^ to be equal to the Borel parameter /i^ = M^. The results take the 
form: 



1 




(22) 








l/2 + x + ln{{l-xf/x) + 



-5 - 17x + 16x2 + 12x^ 



(3-2x)xHn(l/x) 
(l-xy(l + 2x) 



(23) 



6{l-x){l + 2x) 




6(1 -a;)(l + 8a;) (1 - x)2(l + Sx) J J 

where uo{x) and do{x) are bare loop contributions, given by (19). 



(24) 
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4. Power corrections 



In this Section the power corrections to the sum rules are calculated. 

The lowest dimension (d=4) power corrections are the corrections due to gluon condensate 
{Q\{as/'n)G'^j^ G'^jy|0)- Some examples of the corresponding diagrams are given in Fig.3. The 
calculations are performed in the fixed point (Fock-Schwinger) gauge. The corresponding 
expression for quark propagator is given in ref.l2, eq.(24). (The error in the coefficient in 
front of the last term was corrected: it should be 1/288 instead of 1/96). In order to be sure 
in the final results the fixed point was chosen in two ways: at the upper and lower left-hand 
vertices of Fig.l diagrams. The results for the sum of diagrams coincide in these two cases, 
as it should be. However, the contributions of separate diagrams are different. Particularly, 
if the fixed point is chosen at the upper left-hand vertex, then all diagrams where the soft 
gluon is emitted from the upper horizontal line (i.e., by active quark) in Fig.l, are zero in the 
lowest twist approximation. However, if the fixed point is at the lower left-hand vertex, the 
diagrams with soft gluon, emitted by active quark are nonzero. [] Therefore, it is generally 
untrue the folklore statement, that in the lowest twist approximation one may neglect in the 
forward scattering amplitude the active quark interaction between absorption and emission 
of virtual photon - this statement is gauge dependent. 

All other technique is the same as in ref.15. The contributions of gluon condensate to 
u and (i-quarks distribution were found to be (the ratios to bare loop contributions are 
presented) : 

u{x){G^) ^ ((a./Tr)^') (11 + 4x - 31a:') _ /^^2 

uo{x) M4 12 a;(l-a;)2(l + 8a;) ^ NP' ^ ' 

d{x)^G'^} _ {{as/7r)G^) 7i^ (1 - 2a:') sp 
do{x) M4 6 x^{l-xy{l + 2xy ^^M^' ^ ' 

(The factor (1 — e~*o/*^^) appears since gluonic condensate also contributes to continuum 
and this contribution should be subtracted). 

Consider now the contributions of d=6 operators. 

We start from studying the diagrams of Fig.4 and Fig. 5, the contribution of which is pro- 
portional to the vacuum mean values of the type {ipa'ip/sDpG^^) , {ipaV oi'isG^^,) , (ip'V'Wip) 
ets. All of them are expressed through (ipip)^ (see [12]). Here it should be necessary to make 
the following remark: since the approach is inapplicable at x ~ 1, then the diagrams of the 
type of Fig.2, the imaginary part of which is proportional to 6{1 —x), should not to be taken 
into account. But the diagrams of Fig.4 arising by evolution of the diagram Fig.2, should 
not be taken into account too. The diagrams which contribute to intermediate region x (and 
not being evolution of the corresponding non-loop ones) are given in Fig. 5 (cf. Fig. 5b in 
[12]). These diagrams are one-loop ones and for this reason it may be expected that their 
contribution will be dominant among contributions of d=6 operators. The contribution of 
Fig. 5 diagram into sum rules for quark distribution in proton appeared to be equal to: 



^May be, this observation is related to the Brodsky et al. [20] stateraent, that structure functions are 
influent by active quark interactions with spectators. 
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Q;,a2(215 - 867x + 172x2 + 288(1 - x)ln2) 1 



Uq{x) 



M6 • 8l7r • 8a;(l - x)3(l + %x) 



E2(so/M^) 



(27) 



do{x) ~ M*^81nx{l-xf{l + 2x) E2{sq/M^ 



(28) 



where 



a = -{2Tif{{)\t^\{)) (29) 

Besides the considered above contribution of the quark condensate square q;s(0|V^V'|0)2, 
there exist contributions of d=6 operators described by two-loop diagrams. Examples of 
these diagrams are shown in Fig.6. The contributions of these diagrams are expressed via 
a1{Q\'4)'il)\Q)Y and (0|5'^/"*''G^j,G'^;^G'^^|0). The diagrams of the first type are suppressed in 
comparison with (27), (28) by agj-K and by numerical factor ~ l/27r, because they are two 
loop diagrams. The diagrams of the second type, proportional to and (0|5'^/"'"^G^^G^;^G^^|0) 
are also suppressed by two-loop factor l/27r, but vacuum expectation value of operator is 
poorly known: there is only instanton estimate [21] 



487r2 



(l/pD(0|(«./7r)G2|0) 



(30) 



where pc is the effective radius of instanton. The magnitude of pc is not well known: various 
estimations result in essenially different values, varying from 1/3 [22] to 1 fm [21], Also, one 
may have doubts [23] if instantons quantitatively describe physical processes at the scale 
~ IGeV'^, where we are working. Moreover, some of these diagrams are infrared divergent 
(for examples see Fig. 6b). 

For these reasons we restrict ourselves to the estimate of (G^) contribution only. Such 
estimate shows, that it may be essential in the domain of small x < 0.2 and is positive in 
both cases - for u- and (i-quarks. 

The final result for valence quark distribution in proton are of the form 



xu{x) ~ 



.2x{l-x)'{l + 8x)E2{^) 



1 



Uo{x) 



U{X){G^)+U{x)^^(^^^^2 



Uo(x) 



+ 



(31) 



xd{x) = 2a:(l - x)\l + 2x)^2(^) 



1 



do{x) 



1 + 



^^^(^,Qg) 

do{x) 



+ 



(32) 
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5. Valence u- and d-quark distributions at intermediate x. 



We arc now in a position to calculate numerically valence u- and d-quark distributions 
using eqs.(31),(32). The following values of the parameters were taken: X% — 2.1GeV^, 
So — 2.3GeV'^ [24], Xqcd — 250MeV. The latter is the effective one-loop value, which gives 
the same ^^(Q^) as many-loop calculations at low = l — bGeV"^ and particularly, ^^(m^) = 

0. 355 found in [23]. The parameter a,(0|^z/'|o)2 was varied from SAQ-'^GeV^ to 2.2.10-^GeV^ , 

1. e, asO^ varied from O.lSGeV^ to O.SiGeV^. (The lower limit is the old value [25], the upper 
limit corresponds to the recent determination of this parameter from r-decay data [26], [27]). 
The important parameter is the value of gluon condensate b — (0](Q;s/7r)G^^|0). We allow 
it to vary from up to 0.012GeV'^. (The recent determination of (0|(Q;s/7r)G^^|0) from the 
analysis of the r-dccay data results in: (OKa./Tr)^^ jO) = 0.009 ± OMTGeV^ [23], the value 
0.012Ge1/^ is the old result [28]). AU results are presented at Ql ^ 5GeV^. 

Let us start with studying of the Borel mass dependence of quark distributions at various 
X. (Fig. 7; the chosen values of the parameters are: — O.SAGeV^, b — O.OOdGeV^ - our 
favourite values - see below.) As is seen from Fig. 7, for -u-quark the stability in is good 
or satisfactory at 0.1 < a; < 0.5, at a; = 0.6 the dependence is rather strong. For o?-quark 
distribution the stability interval is narrower; 0.25 ^ x < 0.55. 

Check now how much are the contributions of power corrections in comparison with d — 
OPE term (bare loop). For it-quark the applicability domain at low x is limited by gluon 
condensate contribution - it comprises about 30% at a; = 0.15 (and = l.lGeV'^). At the 
same x = 0.15 the contribution of (G^) term becomes large, although uncertain. At large 
X the limits for w-quark come from perturbative corrections, they comprise about 40% at 
X = 0.65. For d-quark the corresponding limits are 0.25 ^ x < 0.60. (All the mentioned 
above values are given for a.a^ = OMGeV^, b = 0.006G'ey^ = l.lGeV^). 

Fig.8a,b shows the dependence of u{x), d{x) on the magnitude of gluon condensate. One 
can see from Fig.8, that the data following from the analysis of deep inelastic scattering are 
best described at 6 = 0.006Gey*^. The zero value of gluon condensate cannot be excluded 
with certainty, but higher values, b > 0.012GeV'^ give much worse description of the data. 
The curves in Fig.8a,b were calculated at a^a^ = 0.34^6^^ (and = l.l.GeV'^).The 
dependence of quark distributions on the magnitude of quark condensate square times - 
a.a^ = (27r)^a,(|V5V^|0)2 is plotted in Fig.9a,b (at b = O.OOGGeV^). The best fit is obtained 
at agO^ = 0.34^61^^. Therefore, as follows from the analysis, our favourite values of gluon 
condensate and times quark condensate square are 

(oi^G^ |o) ^o.ooecey^ 

TT 

a,(0|^^/'|0)2 = 2.2. ■ lO-^GeV"^ (33) 

which agree with the recent result from the r-decay analysis [23], [26], [27]. 

Besides, the uncertainties arising from the spread of possible values of b and a^, there 
are additional sourses of errors - the uncertainty in and those coming from continuum 
contribution. The value of = 2.1GeV^ was found in [24] by the best fit of the proton 
QCD sum rules. However, in [24] the old value of as{0\ipip\0)'^ = 8.10''^GeV^ was used and 
as corrections were not accounted. This may result in 10 — 20% deviation of A^ from the 
above accepted value. Continuum contribution comprises about 30% of the total in w-quark 
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case and about 60% in d-quark case. So, one may estimate possible errors from uncertainty 
of continuum model as 10% and 20%, correspondingly, for u- and d- quarks. Therefore, 
our estimates of total errors are: for li-quarks distributions, in the middle of x-interval - 
0.25 < a; < 0.45, about 25%, at the ends of intervals - x = 0.15 and 0.65 - about 50%; for d 
quarks distributions, in the middle of x-intcrval - 0.3 < a; < 0.45 about 30%, at the ends of 
intervals - x = 0.25 and 0.55 by a factor of 2. Taking into account these errors, the agreement 
with u- and d-quark distributions, found in [4] from experimental data (solid curve in Fig. 8) 
is satisfactory. It should be mentioned, that the account of (G^) terms would improve the 
agreement in the domain of low x. 

A few additional remarks are in order. In Figs. 8, 9 wc performed the comparison with 
valence quark distributions found from hard processes, we used only one group results [4]. 
We hmited our comparison to these data not because we consider these results as more 
confident or rehable, than those of other groups [5-9]. Quite opposite, we believe, that the 
precision of new analysis (sec especially [29], [30]) is better than the older ones. The reason is 
that theoretical errors in determination of u- and d- quark distribution exceed the differences 
in various treatment of the data - all of them are in the limit of our theoretical uncertainties. 
We chosed for comparison the LO results of [4], because our perturbative calculations are 
done also in LO. Our results can be compared directly with experimental data by considering 
the difference 

F^x) - F^{x) - ^[u{x) - d{x)U + l[u{x) - d{x)] (34) 

(x) — -F^(x) was measured by NMC [31]. At a; > 0.30 one may expect, that sea quark 
contribution - the last term in (36) - is negligible and the data can be compared with valence 
quark distributions. The result is that at 0.3 < x < 0.55 the theoretical curve (at b = 0.006 
GeV®) is about 30-40% higher than the data points. Such disagreement is in the limit of 
estimated errors. (At b — the agreement is better.) 

6. Conclusion 

The distributions of valence u- and d-quarks at low and intermediate x were theo- 
retically calculated basing on the first principles of QCD: perturbation theory and operator 
product expansion (OPE). No experimental data and no fitting parameters were used. New 
technique, suggested in ref.l5, was exploited: the double Borel transformation in virtualities 
of two hadronic currents in 4-point function eq.(l), which allows one to kill the background 
non-diagonal transition amplitudes. In OPE the operators of dimension 4 (gluon condensate) 
and dimension 6 (quark condensate square times a^) were accounted. The theoretical analysis 
of the obtained valence quark distributions at = 5GeV^ showed, that w-quark distribution 
is reliable at 0.15 < a; < 0.65, its accuracy is about 25% in the middle of this interval and de- 
creases to 50% at the ends of interval; d-quark distribution is rehable at 0.25 < x < 0.55 with 
an accuracy of about 30% in the middle and is given by a factor of 2 at the ends of interval. 
In the limit of these accuracies the theoretically calculated valence quark distributions are 
in agreement with those found from deep inelastic scattering and other hard processes data. 
The account of contribution, probably, improves the agreement, especially at small x. 
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Figure 1: Bare loop diagrams, corresponding to unit operator contribution for u- and (i-quarks 
(respectively, a) and b)). 
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ure 2: Examples of non-loop diagrams {d — 6) for d- and xi-quarks (respectively a),b)). 
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Figure 3: Examples of diagrams for (i = 4 operator contribution, wavy lines correspond to 
gluon. 
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Figure 4: Diagrams for d = 6 contribution, which are treated as perturbative correction to 
non-loop ones and should be omitted. Other notations as in Fig. 3. (For diagrams b), d), e), f) 
symmetrical diagrams are not shown). 
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Figure 5: Examples of o? = 6 diagrams, which are taken into consideration, non-labelled 
quark propagator means that it can be both of u, d-quarks, other notations - as in Fig. 3. 
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Figure 6: Examples of two-loop diagrams of {d = 6),(a,b). Diagrams b) are examples 
infrared divergent diagrams. 
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Figure 7: Borel mass dependence of sum rules for u- and d-quark distributions at various x. 
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Figure 8: d- and w-quark distribution at various values of gluon condensate {b — 0.012, 0.06 
and GeV^, respectively dotted, solid and dashed lines). Thick solid line corresponds to the 
results of [4]. 
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Figure 9: d- and w-quark distribution at various values of quark condensate a^a^ = 
0.13, 0.23, 0.34 GeV^ in comparison with the [4] result (the curve denoted as GRV). 
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